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Motivated by recent observations of d symmetry breaking in strongly correlated two-dimensional 
electron systems on a square lattice, we analyze this phenomenon within an extended Fermi liquid 
approach. It is found that the symmetry violation is triggered by a continuous topological phase 
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mean-field theory, the structure of a part of the single-particle spectrum violating d symmetry is 
found to be highly anisotropic, with a peak located in the vicinity of saddle points. 
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I. INTRODUCTION 

Experimental studies of strongly correlated two- 
dimensional (2D) electron systems have revealed viola- 
tions of the fundamental symmetries of time reversal and 
C4 rotational invariance inherent in the ground states 
of these systems on a tetragonal latticeji"— Considerable 
theoretical effort has been aimed at understanding the 
nature of these phenomena and identifying their under- 
lying mechanisms. Kivelson, Fradkin, and Emerjsi were 
the first to discuss the case of nematic phase transitions, 
well before relevant experimental data became available. 
Somewhat later, Yamase and Kohnc^ (within the t — J 
model) and Halboth and Metzner^ (within the Hubbard 
model) attributed the breaking of four-fold symmetry 
to violation of a Pomeranchuk stability conditioniS, as- 
sociated with antiferromagnetic fluctuations. An analo- 
gous result was obtained by Valenzuela and Vozmediano 
within an extended Hubbard modelf^ 

As a rule, calculations on the ordered side of the 
implicated second-order phase transition are carried 
out within the mean-field (MF) approachi^^rJi An ef- 
fective Hamiltonian containing a separable interaction 
(i2(p)d2(pi) with order parameter d2{pxTPy) = cospx — 
cos Py is adopted to analyze the onset of C4 symmetry vi- 
olation and properties of phases arising beyond the crit- 
ical point. (Momentum components Px,Py are measured 
in units of the inverse lattice constant.) This approach 
has the advantages of transparency and analytical acces- 
sibility. However, it has noteworthy shortcomings, in- 
cluding sacrifice of translational invariance of the inter- 
action. Furthermore, the structure of the relevant order 
parameter is always postulated in the MF theory; as a 
rule the simplest assumption is made consistent with the 
type of symmetry breaking being considered. However, 
the structure of the new ground state often turns out to 
be quite intricate, such that it cannot be properly de- 
scribed in terms of any single order parameter. As will 
be seen, it is just this situation that emerges in dealing 



with the C4 symmetry violation in question. 

In the scenario proposed here, the system is consid- 
ered to be on the disordered side of an antiferromagnetic 
phase transition; hence the corresponding Pomeranchuk 
stability condition is not violated. With the system sit- 
uated far from the transition point, the fluctuation ex- 
change is readily analyzed and is too weak to gap out 
the single-particle spectrum. It will be shown, however, 
that even if the antiferromagnetic fluctuations are weak, 
their momentum dependence is able to promote a topo- 
logical phase transition associated with disruption of C4 
rotational invariance. 

In Sec. II we adopt the Landau-Migdal quasiparticle 
approach to investigate C4 symmetry breaking in a 2D 
electronic system on a square lattice. A simple model 
with an infinite-range interaction function is employed in 
Sec. Ill to analyze a quasiparticle rearrangement due to 
antiferromagnetic fluctuations. In Sec. IV we present and 
discuss results of numerical calculations for a more realis- 
tic model having a flnite-range interaction. Sec. V is de- 
voted to explanation, within the infinite-range model, of 
the arc structure of the Fermi line observed in many high- 
Tc materials. Our flndings are summarized in Sec. VI. 



II. C4 SYMMETRY BREAKING WITHIN THE 
FERMI LIQUID APPROACH 

Adopting the Landau-Migdal quasiparticle picture, in 
which the physical many-fermion system is viewed as a 
system of interacting quasiparticles, the genesis of C4 
symmetry breaking can be investigated based on the fun- 
damental relatio n -'^^d^ 

where dv — dpxdpyj {I'k)'^ is an element of 2D momen- 
tum space. This relation connects the quasiparticle spec- 
trum e(p) with the quasiparticle momentum distribu- 
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tion n(p) = [1 + exp ((e(p) — /i)/T)]~"^ through a phe- 
nomenological interaction function T . This function, 
which is defined by a specific static hmit of the quasi- 
particle scattering amphtude with initial and final ener- 
gies on the Fermi surface, --'^^ depends only on the mo- 
menta p, pi of the colliding quasiparticles. Of the two 
particle-hole channels relevant to the scattering ampli- 
tude T , denoted t and u in the Mandelstam's terminol- 
ogy, the transverse t channel carries vital information in 
the momentum transfer q = P — Pi, whereas the lon- 
gitudinal u channel is silent because the corresponding 
momentum transfer is close to zero. 

In homogeneous matter where total momentum is con- 
served, the first term on the right side of Eq. (H]) is just 
the bare velocity p/M, with M the free particle mass.^° 
In the presence of a crystal-lattice field, the bare group 
velocity is multiplied by a quasiparticle effective charge 
eq(p). However, this modification will be ignored, since 
it reduces merely to a renormalization of phenomenolog- 
ical coefficients ti specifying 2D tight-binding electron 
spectra 
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We are concerned specifically with the impact of an- 
tiferromagnetic fiuctuations on the electron spectra e(p) 
calculated using Eq. ([!]) . Treatment of the effect of these 
fluctuations on the interaction J- does not encounter dif- 
ficulties far from the attendant antiferromagnetic phase 
transition. The corresponding fluctuation exchange is 
adequately addressed within the Ornstein-Zernike (OZ) 
approximation, which neglects scattering of fluctuations. 
The part of J- responsible for the exchange is then 

Kp-ysiP, Pi) = X^CTapCr^S [(p - Pi - Q)^ + C^] • 

(3) 

The constant A represents the spin-fluctuation vertex and 
Q = (tt, tt) is the antiferromagnetic wave vector, while ^ 
is the correlation radius. 

Inserting Eq. ^ into Eq. ((!]) and evaluating the spin- 
fluctuation contribution aided by the identity 2(Tai3(Tjs = 
SSasS-yp — (TaS(^-f0, One arrives at 
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The normalization condition 2 J n{p)dv — p determines 
the chemical potential /i consistent with density p. This 
approach to the problem is self-consistent provided the 
dimensionless parameter fN{0) is rather small, where 
/ = (3A^/47r) ln(l/^) is a coupling constant and N{0) ~ 
l/27rio is the density of states of a 2D electron gas on a 
square lattice having the tight-binding spectrum ([2|). 

Direct numerical solution of this 2D nonlinear integral 
equation is extremely time-consuming. If only the com- 
ponent of the interaction ([3]) proportional to (i2(p)d2(Pi) 
is retained, then beyond the point where the correspond- 
ing Pomeranchuk stability condition is violated, one ob- 
tains the ordinary mean-fleld theory equations. However, 



this approximation is quite poor for the interaction ([3]), 
which peaks at momentum transfer q = Q. Accordingly, 
the customary MF scenario must be regarded as vulner- 
able. 

Our approach to the problem stems from this observa- 
tion: collapse of collective degrees of freedom associated 
with violation of sufficient conditions^" for the stability 
of the standard Landau Fermi Liquid (FL) state is not 
the only possible scenario for the breakdown of C4 sym- 
metry. A viable alternative is provided by violation of 
a necessary stability condition. This condition requires 
that an arbitrary admissible variation dn{p) from the FL 
quasiparticle momentum distribution np{p), while con- 
serving particle number, must produce a positive change 
of the ground-state energy Eq, 
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(e(p; riFip)) - p) 5n{p)dv > 0, 



(5) 



where e(j);np) is the spectrum of single-particle excita- 
tions and p the chemical potential. 

Violation of the condition ^ is accompanied by a 
change of the number of roots of the equation 



e{p,nF) = p, 



(6) 



which implies a change of the topology of the Fermi sur- 
face. For a thorough development of the concept, see 
the review by Volovik,^- Throughout, we adhere to his 
rigorous quantitative definition of topological phase tran- 
sitions, as distinguished from looser notions such as tran- 
sitions between large and small Fermi surfaces that are 
also prevalent in the literature. It should be emphasized 
that in contrast to the original Lifshitz description,— the 
topological transition under consideration is triggered by 
the interaction between quasiparticles (see also Refs. |2J- 

m. 



III. QUASIPARTICLE REARRANGEMENT 
WITHIN A SIMPLIFIED MODEL 



To gain insight into the essence of this scenario, we 
restrict the analysis to zero temperature and simplify the 
interaction. Replacement of the interaction term (j3|) by 
an infinite-range form ~ i5(q — Q) leads directly to the 
explicit versionSi 



e(p)=e° +/n(e(p + Q)) 



(7) 



of relation ([T]), where / is the coupling constant identi- 
fied above. This treatment is analogous to that adopted 
by Nozieres^ in a study of non-FL behavior of strongly 
correlated Fermi systems in the case where forward scat- 
tering in the t channel prevails. Eq. ^ can be derived 
within a standard variational procedure based on the 
formula^! 



J [eln{p) + i/n(p)7i(p + Q)]2dv 



(8) 
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for the energy E of the model quasiparticle system. This 
form for the energy functional admits a greatly simplified 
analysis of the problem due to the partial separation of 
different p channels. 

To proceed, we observe first of all that at T = 0, the 
posed rearrangement of the initial standard Landau state 
can occur only in those 2D systems where there exist hot 
spot&2? — points situated on the Fermi line and connected 
by the vector Q. Indeed, in systems with small quasipar- 
ticle filling, the product n(p)n(p + Q) vanishes for any 
momentum p, so that the ground-state energy is inde- 
pendent of the coupling constant /. The same is true in 
the case of small hole filling. 

In systems with hot spots, the rearrangement oc- 
curs due to breaking of the quasiparticle pairs occupy- 
ing single- particle states with momenta p and p -|- Q. 
The corresponding domain TZ (the "reservoir" ) consists of 
four quasi-rectangles, each adjacent to one of the saddle 
points (0, tt), (tt, 0), (0, — tt), (— tt, 0) of the tight-binding 
spectrum 6°. Each of the four elements of TZ is con- 
fined between (i) the border of the Brillouin zone, (ii) 
the counterpart of the initial Fermi line, defined by the 
equation Cp+q — fJ-, and (iii) two segments of the Fermi 
line embracing the given saddle point. 

Quasiparticles move out the domain TZ to resettle in 
a region C where all pairs of single-particle states con- 
nected by the vector Q are empty. The region C com- 
prises four "lenses," situated between neighboring hot 
spots and bounded by the initial Fermi line and its coun- 
terpart (see panel (a) of Fig. [T]). The transfer of one 
quasiparticle from TZ to C produces a gain in energy 
which is just the coupling constant / minus the loss r 
of kinetic energy. Its minimum T,„in is attained when 
a quasiparticle, vacating a state in TZ with momentum 
p, occupies in £ a state of lowest energy, given by the 
chemical potential, so that r^in = ~ £p- Therefore the 
rearrangement is favorable provided — ^ + f > 0. 

An alternative process involves transfer of the quasi- 
particle counterpart, which has momentum p+Q. In this 
case, the rearrangement occurs provided Cp+q— > 0. 
The choice between the two options is decided by compar- 
ing the corresponding energies. The boundary at which 
one behavior gives way to the other is determined by the 
relation — Cp+q- Since the straight line so defined is 
part of the new Fermi line, we infer that the rearrange- 
ment has converted the original, isolated hot spot into a 
continuous straight line of hot spots, i.e., a hot line (HL) 
(see panel (b) Fig. [T]). 

These results imply that quasiparticles are swept from 
a certain subdomain 5 of 7^ consisting of eight approx- 
imately trapezoidal strips. The boundaries of a given 
strip are traced on three sides by (respectively) the ini- 
tial Fermi line, the border of the Brillouin zone, and a 
line geometrically similar to the initial Fermi line but 
shifted into the domain TZ (see Fig.[T]). The strip's fourth 
side (red on-line) is just the hot line. This solution is 
self-consistent: any single-particle state with momentum 
p e iS has its counterpart, with momentum p+Q, located 



outside S, and this state is occupied, so that Eq. ([7]) is 
fulfilled. Transparently, the new momentum distribution 
does not violate C4 symmetry. 

Defining the strip energy width Ws of the region S as 
the maximum of the initial hole energy [Cp — /i^j consis- 
tent with the rearrangement, one has Ws — Di ~ Df, 
where 2Di (respectively, 21?/) is the minimum energy 
distance between the segments of the initial (final) Fermi 
line situated in different half-planes. On the other hand, 
one finds Ws — f — {fi — fj-i), where fii is the initial 
chemical potential. To estimate the strip width Ws and 
the difference 11 — fii, both proportional to /, we (i) ap- 
proximate the Fermi velocity v''(p) — (^d€p/dp)^ on the 
Fermi line in terms of two parameters, namely its av- 
erage magnitudes vf and in the lens and strip re- 
gions, respectively, and (ii) invoke the coincidence of the 
chemical potential with the Fermi energy that is intrin- 
sic to Landau theory. In the lens region C one then has 
fi — fii — vfwi, where wi is the momentum width of the 
lens filling. In the domain S, one obtains the analogous 
formula Ws = f — (l^ — IJ^i) = WgVs- Particle-number con- 
servation implies that wik = 2ws{ls — Ws/i), where Is is 
the strip length, li is the lens length, and Ws = Ws/v^s 
the momentum width of the strip. Upon elimination of 
Ws and wi from these relations, we arrive finally at 
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(9) 



for small Ws- 

As long as all the saddle points remain occupied, C4 
symmetry is preserved. However, as the electron density 
p decreases, the distance between the new Fermi line and 
the saddle points shrinks. At a critical density pc, or 
equivalently, at the critical constant fc where two seg- 
ments of the Fermi line that cross the same boundary of 
the Brillouin zone merge at the saddle point, the number 
of solutions of Eq. ([6]) certainly drops, thereby signaling a 
topological phase transition. In the critical situation one 
has Df = 0, or equivalently Ws = Di. The trapezoidal 
shape <S then becomes triangular, and we have 
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Using this result, the critical value Fc of the dimension- 
less constant F — fN{0) is given by Fc = D/2ntQ. As- 
suming the ratio D/to to be small, we thus have Fc ■^l, 
which implies that the derivative i9E(p,e)/9£ remains 
small, i.e., that the e-dependence of the mass operator 
S(p,e) is moderate.'^"^ Under these conditions, the gen- 
eration of new branches of the single-particle spectrum 
e(p), such as the small pockets of the Fermi surface sug- 
gested to explain magnetic oscillations in the pseudogap 
regime?^ is questionable. 

Beyond the transition point (e.g. at p < pc), C4 sym- 
metry is necessarily broken. Suppose, conversely, that it 
is preserved. Then all the saddle points must then be 
emptied simultaneously, implying that every rearranged 
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FIG. 1: Color online. Panel (a): Fermi line (black) and 
its counterpart (blue(gray)) for the bare tight-binding spec- 
trum Ep = —2to (cospx + cospy) ~\- 4ti cospxcospy, with 
ti/to = 0.45. The reservoirs TZ are colored in black, and 
the lenses £., in light gray. The hot spots connected with 
each other by the vector Q are symbolized by red (gray) 
dots. Panel (b): Fermi line for the model assuming the 
infinite-range interaction function /(q) = (27r)^/(5(q — Q), 
with fN{0) — 0.13. Hot lines are drawn in red (gray). Fermi 
lines for the bare tight-binding spectrum and its counter- 
part are shown as green (light gray) and blue (gray) lines, 
respectively. 

saddle point energy Cs [p) exceeds the chemical potential 
/u(p). But according to Eq. ([7]), the interaction con- 
tribution to ts vanishes when all the saddle points are 
emptied. Hence the saddle-point energy es{p) must co- 
incide with the corresponding bare value e^{p), implying 
that egip) > p{p)- However, if the difference pc — P is 
small, then without fail e1{p) < Pi{p). Thus a contra- 
diction is encountered, since it follows from Eq. (|9]) that 
Pi{p) < pip)- This deadlock is resolved if, beyond the 
critical point, only one of two neighboring saddle points 
is emptied, with the second remaining occupied. Such a 
solution is indeed consistent with Eq. ([7]). 



IV. NUMERICAL RESULTS WITH REALISTIC 
INTERACTION 

The results we have derived for the simple model based 
on an infinite range interaction ~ (5(q — Q) are in agree- 
ment with those obtained from numerical calculations 
performed for the more realistic interaction (jS]) and dis- 
played in Figs. [2] and ([3]). Some complications associ- 
ated with the finite correlation radius of the interaction 



^ will be considered below, but first we examine the 
results of the extended Fermi-liquid theory in compari- 
son with corresponding predictions of MF theory. The 
MF single-particle spectrum coincides with a bare spec- 
trum before the transition point is reached, while be- 
yond the transition it receives a correction JemfIp) — 
77 (cos — cos py)^ with the order parameter rj taking the 
same value throughout the Brillouin zone. The Fermi 
line calculated within the extended FL approach devi- 
ates substantially from that predicted by MF theory. In 
particular, upon comparing the upper and lower panels 
of Fig. [21 we see that in the lens domain the location 
of the Fermi line remains almost unchanged as the sys- 
tem passes through the transition point. Indeed, this 
behavior also prevails over a significant portion of the 
HL region away from the saddle points. In other words, 
beyond the point where C4 symmetry is lost, the associ- 
ated rearrangement of the Fermi surface occurs only in 
the immediate vicinity of the saddle points — in a sharp 
contrast to what is found in MF theory. 

Analogous conclusions follow from a study of Fig. |31 
where the Fermi velocity calculated on the basis of Eq. ([T]) 
is plotted. It is seen that the correction to the the bare 
Fermi velocity Vp stemming from antiferromagnetic cor- 
relations as described by ^ remains smooth and small 
except in the HL region, where it soars upward. 

Such behavior of the Fermi velocity vp, which per- 
sists through the transition point, can be elucidated by 
analyzing the Landau relation ([T]). First, we observe 
that the overwhelming contributions to the integral in 
this relation come from the HL region; otherwise there is 
no appreciable overlap between the peak in the interac- 
tion function and the 6 peak in the derivative dn{p)/dp. 
To proceed further, we introduce a new set of orthog- 
onal momentum coordinates pt,Pn, with the axis pt di- 
rected along the HL and the axis Pn perpendicular to 
it. In the HL region we then have dp^dpy = dpndpt 
and dn{p)/dpn = ±(5(p„), the sign of the derivative be- 
ing positive in the left half-plane and negative otherwise. 
This alternate sign is responsible for the vanishing of the 
group velocity at the saddle points, through interference 
of the contributions to the integral term in Eq. ([IJ from 
neighboring segments of the Fermi line situated in the two 
half-planes. The distance between these segments (as de- 
fined in Sec. HI) is 2Df. If the inverse correlation radius 
which measures the radius of the spin-interaction 
term ([3]) in momentum space, turns out to be so small 
that ^"^ < Df, then the two contributions cease to in- 
terfere, and the elevation of the HL value of the Fermi 
velocity is readily estimated as 

VFipt) ~ (27r)-2r^ (11) 

These conclusions are in agreement with the results 
for the model with finite-range interaction presented in 
Fig. [3l At the same time, the estimate ([TT|) is in agree- 
ment with the jump of the single-particle spectrum e(p) 
on crossing the HL, found for the simple model with 
(5(q — Q) interaction and implying an infinite value of 



the model's HL group velocity. The above considerations 
demonstrate that the FL rearrangement of the ground 
state leading to the phenomenon of C4 symmetry viola- 
tion has little in common with the rearrangement pre- 
dicted by conventional MF theory based on the single 
order parameter ^2- 

The analysis can be made more informative by focusing 
on the difference D(px,Py) = t{px,Py) — <^{Py,Px) and its 
integral D over the intermediate momenta Px,Py Both 
quantities vanish on the disordered side of the phase tran- 
sition, and beyond the transition point it is straightfor- 
ward to evaluate D by means of Eq. Q. For — > 0, 
one can make use of the formula n{px,Py) — n{py,Px) — 
{dn{p/d€{p)D{px,Py) to recast this equation in a form 

equivalent to the Pomeranchuk stability condition, whose 
violation is a prerequisite for the MF description of C4 
symmetry breaking. From the preceding discussion, we 
infer that if a nontrivial solution of Eq. ([T^ exists, it 
must be anisotropic, with a peak located in the HL do- 
main and having a width of order Such a structure 
of the order-parameter function D(j)) is quite unlike that 
adopted in conventional MF theory of the observed x — y 
symmetry violation. In evaluating the integral in Eq. 
we employ the relation dn{p)/de{p) = {dn{pn)/dpn)/vp. 
Referring to the above derivation of the estimate ((TT|) , it 
is seen that the relevant value of the group velocity is 
vp — ^, as long as ^ > {Di — Ws)~^. The ^ dependence 
of the integral is then effectively nullified, precluding non- 
trivial solutions of Eq. (|12p . 

Nontrivial solutions of Eq. (|12p can in fact emerge be- 
fore the two neighboring pieces of the Fermi line meet 
each other at the saddle point, provided ^ < {Di — Ws)~^- 
In this case, the characteristic value of the Fermi ve- 
locity drops somewhat, thereby enhancing the integral. 
Whether this enhancement is sufficient for the violation 
of the Pomeranchuk stability condition will be decided in 
a more intensive round of numerical calculations. 

It is worth emphasizing that the situation underly- 
ing the violation of C4 symmetry in systems in which 
the Fermi surface comes close to van Hove points is not 
specific to either the MF treatment or our analysis. In 
fact, the effective Stoner factor, which determines the 
enhancement of the effective field acting on a particle in 
matter, is proportional to the product of the interaction 
strength and the density of states. The latter diverges 
at a van Hove point, and hence the corresponding Stoner 
factor diverges as well, independently of the shape of the 
order parameter. The crucial point of distinction is as 
follows. In MF theory, which reasonably exploits the en- 
hancement of the density of states near the van Hove 
points and an order parameter d2{p) having the needed 
symmetry, the effective field stretches over the whole Bril- 
louin zone in accordance with the chosen shape of the 
order parameter. In our approach based on exchange of 
antiferromagnetic fiuctuations between electrons, it is in- 
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FIG. 2: Color online. Fermi lines for the model assuming the 
finite-range interaction function /(q)=/a/((q— Q)'^-l-C~^), 
with ^=30. Panel (a): /aA^(0)=0.32; C4 symmetry is not bro- 
ken. Panel (b): /£,A''(0)=0.48; one of the two solutions with 
spontaneously broken C4 symmetry is shown. Only the first 
quadrant of the Brillouin zone is drawn since neither p^^—Px 
nor Py— > — Py reflection symmetry is broken. Fermi lines for 
the bare tight-binding spectrum ep and its counterpart are 
shown as green (light gray) and blue (gray) lines respectively. 



stead the shape of the exchange interaction that governs 
the behavior of the effective field. Since this field dies 
out at rather small distances from the saddle points, the 
topological rearrangement of the Fermi surface violating 
C4 symmetry occurs only in the regions close to these 
points. 

Let us now identify inherent properties of the interac- 
tion function J- responsible for violation of C4 symmetry 
and more generally for topological transitions. In ho- 
mogeneous matter, it is well understood that topological 
phase transitions are characterized by a change in the 
number of sheets of the Fermi surface.— In 2D electron 
systems on a square lattice, topological transitions are of 
much the same character. The salient common feature 
here is that no symmetry is violated, provided that a local 
rearrangement of the quasiparticle momentum distribu- 
tion leads to dominance of forward scattering in the t 
channel referred to the momentum transfer q specifying 
J^(q). On the other hand, in the case of antiferromag- 
netic fiuctuations backward scattering prevails. Then, at 
the transition point, quasiparticles leaving the vicinity of 
one saddle point may move into the vicinity of a neigh- 
boring saddle point. Thus the sheet number remains un- 
changed; instead, the symmetry of the ground state is 
violated. 
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FIG. 3: Color online. Fermi-velocity magnitudes vp = 
\de{p) /dp\ (in units of 2to), evaluated along the Fermi line 
as a function of the angle Lp defined in the inset, for differ- 
ent single-particle spectra e(p). Results are shown for the 
bare tight-binding model with the same parameter choice as 
in Fig. [T] (brown (dark gray) line) and for the Fermi-liquid- 
theory model of Fig. [2] at faN{0) = 0.32, T = 10"" (green 
(light gray) line); faN{0) = 0.48, T = 10""* (red (black) line); 
and /aA*'(0) = 0.48, T = 10"^ (blue (gray) line). Broken C4- 
symmetry of the solid/red curve with respect to x—y exchange 
is manifested by its different behavior in the two shaded areas 
close to the saddle points. 



V. ARC STRUCTURE OF THE FERMI LINE 

The model we have developed may also have a bearing 
on the emergence of the arc structure of the Fermi line ob- 
served in many high- Tc materials. If we consider pairing 
based on the interaction ([3]), then Eq. ([SJ must be sup- 
plemented by a pairing term (//2)k(p)k(p + Q),^® where 
K{p)={a~^{p)a~^{—p)) is a superfluid density. With this 
modification, Eq. ([7]) as written remains unchanged, but 
the quasiparticle occupation number n(p) acquires the 
BCS form 7i(p)=l/2— e(p)/2i?(p), with quasiparticle en- 

1 /2 

ergy E{p)= [e^(p)-|-A^(p)] . The additional equation 

determines the gap function A(p). In advance of the 
topological phase transition, where C4 symmetry is pre- 
served, a standard nonzero solution of Eq. has the 
property A(p) = — A(p-fQ) exhibited by D pairing, and 
we find 

E{p)E{p + Q) ^ p 

tanh(£;(p)/2T)tanh(£;(p + Q)/2T) 4' ^ ' 

As seen from Eq. (|14p . the associated gap -Emin in the 



single-particle spectrum is suppressed near the diagonals 
of the Brillouin zone, where 

i?n.i„(r = 0)^i^, (15) 

Wi being the total energy lens width. On moving along 
the Fermi line toward the hot line where one has £'(p) — 
E{p + Q), the gap soars upward, with Eq. (fT4|) yielding 

E{p,T = 0)c^^. (16) 

It is important to note that in the HL region itself, the 
gap value is markedly suppressed, because Eq. (O tells 
us that |e(p)| ~ / in a significant part of this region, 
which is incompatible with Eq. (jl6|) . This indicates that 
pairing has little impact on the violation of C4 symmetry, 
which primarily involves the immediate vicinities of the 
hot lines. 



VI. CONCLUSION 

In summary, we have addressed the problem of C4- 
symmetry violation in electron systems on a square lat- 
tice within a self-consistent Fermi liquid approach, as- 
suming that the Landau interaction describes the ex- 
change of antiferromagnetic fluctuations, which is treated 
within the Ornstein-Zernike approximation. We have 
demonstrated that as the strength of this interaction 
builds up, the distance between saddle points and the 
Fermi line shrinks, eventually generating a quantum criti- 
cal point of a new type, at which a continuous topological 
phase transition triggers the violation of C4 symmetry. 
The group velocity becomes finite again once the transi- 
tion point is passed. Thus, the properties of the electron 
system are governed by Fermi-liquid theory throughout 
the vicinity of the proposed quantum critical point, im- 
plying that magnetic oscillations should be observed on 
both the sides of the topological transition, in agreement 
with recent measurements 1^ 
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